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Abstract 

We consider a one-parameter family of matrix product states of spin one 
particles on a periodic chain and study in detail the entanglement properties of 
such a state. In particular we calculate exactly the entanglement of one site with 
the rest of the chain, and the entanglement of two distant sites with each other 
and show that the derivative of both these properties diverge when the parameter 
g of the states passes through a critical point. Such a point can be called a point of 
quantum phase transition, since at this point, the character of the matrix product 
state which is the ground state of a Hamiltonian, changes discontinuously. We 
also study the finite size effects and show how the entanglement depends on the 
size of the chain. This later part is relevant to the field of quantum computation 
where the problem of initial state preparation in finite arrays of qubits or qutrits is 
important. It is also shown that entanglement of two sites have scailing behavior 
near the critical point. 
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1 Introduction 



Interacting spin systems are among the most promising candidates for the actual 
implementation of quantum computers in the future. In the new terminology which 
has emerged since the upsurge of interest in quantum computation, a spin 1/2 system 
refers not only to the actual spin degrees of freedom of an atom or nucleus, but it refers 
to any system, with any number of levels, in which we have selected two states for 
encoding the information. For example the ground state and the first excited state of 
an ion in an ion trap, make a spin 1/2 system or a so called qubit. Similar terminology 
is used for three level states or qutrits. These latter systems and their generalizations 
to d— level states or qudits are of immense interest, since it is not yet clear if actual 
implementation of quantum computers will be based on two level systems. In view 
of this, many systems which have been traditionally the focus of study in condensed 
matter physics are being examined from different points of view related to quantum 
computation and information. The most important property of an interacting spin 
system, which is of relevance to this new emerging field, is entanglement or non-local 
quantum correlation. In fact it is considered as a resource, like energy, since it plays 
a vital role in any process of quantum information and computation, moreover it can 
be measured, manipulated and transferred. Consider the ground state of an inter- 
acting spin system, comprising spins. There are some basic questions regarding 
entanglement: How much two distant spins are entangled with each other? How this 
entanglement varies with the system size A'^? Is there any threshold distance, beyond 
which there will be no entanglement at all? How the entanglement varies when we 
approach a point of quantum phase transition? Answering these questions requires 
tools which have been developed only recently in the field of quantum information 
[D El El m El E] . In this way a fruitful field of investigation at the borderline of con- 
densed matter physics and quantum information has emerged. The aim of this letter 
is to investigate such questions for a class of spin-one states, which are known to be 
exact ground states of certain multi-parameter families of Hamiltonians describing 
nearest neighbor interactions of spin-one particles on a periodic chain. The method 
for construction of such states, known also as the matrix product formalism, was first 
introduced in [71 El E] and then applied to various models in [lOl HH [HI [13]. Re- 
cently it has been applied even to two dimensional models [HI [15]. Such states can 
be constructed to have specified symmetry properties or even to induce certain kinds 
of quantum phase transitions with pre-determined properties [16]. The entanglement 
properties of the so called AKLT models [7], which inspires the matrix product states 
was first studied in [17j. Here we study the entanglement properties of a one param- 
eter deformation of AKLT models. 

We first study the entanglement of one site with the rest of the chain and the entan- 
glement of two distant sites in the thermodynamic limits, however we mainly focus on 
states with finite but arbitrary number of spins, since this is the case of interest from 
the point of view of quantum computation and information. We will determine the 
entanglement of two spins in the lattice as a function of their distance, the coupling 
constant of the state or the Hamiltonian, which we denote by g, and the system size. 
The results are that: 1- When the parameter g approaches its critical value go = 0, 
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the range of entanglement increases indefinitely at the cost of its magnitude, fig. ([T]), 
2- for any non-zero value of g, there is a threshold distance beyond which there is no 
entanglement between spins, fig. ([3]), 3- for any two spins with a fixed distance, there 
is a threshold system size, above which entanglement vanishes, fig. and there is a 
scaling behaviour in the entanglement of two adjutant spins which is shown in figure 

dZD- 

Let us first remind the matrix product formalism in a language which we find conve- 
nient. 



2 Matrix Product States 

For a homogeneous ring of N sites, where each site describes a d— level state. The 
Hilbert space of each site is spanned by the basis vectors i = 0, - ■ ■ d—1. A state 

l*> = V'iii2-ijvKl>^2,---,«Af> (1) 

is called a matrix product state, if there exist matrices Ai, i = - ■ ■ d—1 (of dimension 
D) such that 



where Z is a normalization constant equal to Z = tr{E^) and E := X]f=o ® 

The correlation functions are readily calculated in this formalism. For example, 
for the one-point functions we have 

where Eq = J2i j {"^10^)^1 ® Aj. In the thermodynamic limit, only the largest eigen- 
value of E survives and so any level crossing in the largest eigenvalue of E, leads to a 
discontinuity of correlation functions. This may be termed an MPS-quantum phase 
transition [T6l. 



2.1 Gauge Transformations 

From ([2]), it is evident that two sets of matrices {Ai} and {A[} lead to the same matrix 
product state if they are related as A'^ = ^SAiS~^, where /i is a scale factor and 
5 is any invertible matrix. Actually the gauge transformation can be more general 
than this, namely A'^ = fiSAiS' with S'S = I. Such transformations can be used to 
gauge away irrelevant parameters in the matrices Ai. 



2.2 Symmetries 

On a ring, the state ([1]) is invariant under translation. Demanding more symme- 
tries imposes constrains on the matrices Ai. Considering equation Q, the state is 
symmetric under parity if there exists a matrix IT such that 

Af = aUAiU-^ a = ±l (4) 
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where A is the transpose of A and it has time reversal symmetry if the matrices Ai 
are real. 

Consider now a local symmetry operator R acting on a site as R\i) = Rji\j) where 
summation convention is being used. R is a d dimensional unitary representation of 
the symmetry. A global symmetry operator TZ := R^^ will then change this state to 
another matrix product state 

M/,,,,...,^^vf':=tr(4^^:^---A^^), (5) 

where A'- := RijAj. The state |^') is invariant under this symmetry if there exists an 
operator U {R) such that 

R^jAj = U{R)AiU-^{R). (6) 

Thus R and U{R) are two unitary representations of the symmetry, respectively of 
dimensions d and D. In case that i? is a continuous symmetry with generators T^, 
equation dH), leads to 

{Ta)^,Aj = [%.A,l (7) 

where Ta and Ta are the d— and D— dimensional representations of the Lie algebra of 
the symmetry. 



2.3 The Hamiltonian 

To construct a Hamiltonian with nearest neighbor interaction, in a way that the state 
in equation ([1]) be its ground state, we have to find the null space of the density 
matrix of two adjacent sites which is given by: 

_tr{{AlA*®A,Ai)E^-^) 
PirM • (8) 

The null space of this reduced density matrix include the solutions of 

J2 CkiAkAi = 0. (9) 

k,l 

The number of independent solutions of this system of equation is d"^ — D^. Thus for 
this density matrix to have a null space it is sufficient that d > D. Let the null space of 
the reduced density matrix be spanned by the orthogonal vectors \ea), (a = 1, • • • s > 
(p _ /)2-j_ Then we can construct the local hamiltonian acting on 2 consecutive sites 
as h := J2q=i ^a\ea){ea\, where A^'s are non-negative constants. These parameters 
together with the parameters of the vectors \ea) inhertited from those of the original 
matrices Ai, determine the total number of coupling constants of the Hamiltonian. 
The full Hamiltonian on the chain is written as H = J2iLi where is the 

embedding of h into sites i and i + 1. The state 1^*) is then a ground state of this 
Hamiltonian with vanishing energy (for an exposition see |15)). 
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3 Matrix Product States on Spin 1 Chains 



3.1 Construction of the state 

The matrix product state, that we use for spin-one systems, has already been con- 
structed in [To]. Here we review its construction in the language introduced in previ- 
ous section for completeness. Since d = 3, to guarantee a null space for the two-site 
density matrix, we set D = 2. The matrices Ai, Aq and Aj = A-i correspond to the 
local states |1), |0) and |1) = | — 1) respectively, where Sz\m) = m\m). Considering 
equation ([7]), the symmetry around the z axis requires that 

[S„Am]=mA,n, (10) 

where 'S'z = ^ ^ q ) Solving ([10]) , demanding parity symmetry ^ and getting 
rid of irrelevant parameters with suitable gauge transformations, leave us with: 

where we have used the compact notation A := X^j^iN) • Note that the state con- 
structed in this way automatically has spin-flip symmetry, i.e. XAmX~^ = a Am with 

X = ^ ^ 0^^ ^ ■ (5)^) = (1> the so-called AKLT point [^, the model 

has full rotational symmetry, since in this case the states — \/2(T+ , cr^ and \/2cr_ , form 
a vector under the adjoint representations of the rotation group. 



3.2 The Hamiltonian 

As mentioned previously, we must solve ([9]) for matrices (jlip to construct the hamil- 
tonian. It is straightforward to verify that the solution space of ([9]) is spanned by the 
following vectors: 

= |1,1) 
les) = -^(|l,0)-a|0,l)) 

1^3) = ! , , (|lJ) + 2g|0,0) + |T,l)) (12) 

|e4> = -^(|0,T)-a|T,0))) 
|e5) = |T,T). 
With these vectors, we write the local Hamiltonian as 

h = a(|ei)(ei| + |e5)(e5|) + 6(|e2)(e2| + |e4)(e4|) + c|e3)(e3|, (13) 
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to preserve the Z2 symmetries mentioned above. Writing this in terms of local spin 
operators, the final form of is obtained as: 

N 

-f^ = ^ ■ Sj+l + J2(Si • Sj+i)^ + Js'S'z.j'S'z.i+l 

i=\ 

+74(5*2,452,4+1)^ + J552,j + JelSj • Sj+1, 52,j52,i+i}+, (14) 

where 

Ji = -6(7(1 + 2g2), J2 = c, 

J3 = (a + H(l + 252), (15) 

J4 = (a + 26(a-l))(l + 2ff2) + (l + 2g)2c, 
J5 = 26(1 + 2^2) + 2c(l-V), 
Je = _6a(l + 2g2)-c(l + 25). 

In writing the above Hamiltonian we have ignored an overall additive constant and 
have re-scaled the operator (jl3p by a constant 2(1 + 1g^^. This represent a four 
parameter family of Hamiltonians which have ([T]) as their ground state. Full rotational 
symmetry exists ,when we have (5,0") = (1, —1), and a = 6 = c = 1, for which case 
the Hamiltonian is known as the AKLT model. 



4 The thermodynamic Hmit 

We can derive many properties of the ground state using the transfer matrix formalism 
explained in section ([2]). The eigenvalues of the matrix E are found to be 

Ai = l + 25 \2 = \-2g A3,4 = fT. (16) 

This shows a level crossing in the largest eigenvalue of E and hence a singularity 
in correlation functions at (7 = 0. The average magnetization and the correlation 
functions are found to be [lOj 

< S\ >=< Si >=< Si >= 0, 

and 

45' fi-M 



< sis: >-- 



(1-2|(7|)2 V1 + 2I5I 
<Sis:,>=-2\g\ia-Sign{g)) 



.1 + 2|5| 

where n is any unit vector in the xy plane. The longitudinal and transverse correlation 
length, diverges at 5 = 0, It is a natural question to ask if the same thing happens 
for entanglement when g approaches this critical point. In [18], the one-site entropy 
which measures the entanglement of one site with the rest of the lattice and also the 
two-site entropy which measures the entanglement of these two sites with the rest of 
the lattice were calculated in the thermodynamic limit. Here we use the negativity 
to measure how much two distant spins are entangled with each other. Moreover, we 
study in detail finite size effects to see how various properties of entanglement depend 
on the system size. 
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4.1 Entanglement of two distant sites 



In the ground state of (|14|) . any two particles will be in a mixed state. The reduced 
density matrix of two spins located at sites 1 and r is denoted by p^'^ . The rotational 
symmetry around the z axis shows that Pij-ki = unless i + j = k + I. Also the 



parity symmetry entails the condition Pij^ki 
calculation shows that 



Pji.-n, where i = —i. Straightforward 



( « 



/3(l,r) 



Hi ^ 

/3 b 
V 6 



\9h 



(17) 



in which 




r-2\ 



0, 7 



a 

' ' VAi 



a J 



1 

Af' 



(18) 



where Ai = 1 + 2\g\ and A2 = 1 — '2\g\. Since the state of two sites, is a mixed state, 
we can not use von Neumann entropy to measure the entanglement between these 
two sites. In |19] it is shown that the necessary condition for a mixed state p, to be 
separable is that its partial transpose has non-negative eigenvalues. The quantitative 
version of this criterion is Negativity which is defined as follows [20] : 



E{p{l,r)) 



\P 



1 



(19) 



where p^^{l^r) is the partial transpose of /o(l,r) with respect to the subsystem A, 
and = ^fX^X is the trace norm of X. Equivalently it is equal to the sum of 

absolute values of negative eigenvalues of the matrix /?^'*(l,r). In other words by 
means of Negativity we can measure the degree to which the partial transpose of the 
state p fails to be positive or how far the state of the two particles is from a separable 
state. A basic property of Negativity is that it is an entanglement monotone, meaning 
that the more entangled a state, the more negative it is in the above sense. 
The eigenvalues of p^^{i,r) are found from (I17p to be: 



^6,7 



"^2,3 = /? 



'^4,5 



\9\1 



\g\l - (5, ^8,9 



a + 7 + 1 



From these eigenvalues the negativity and hence the entanglement of two spins can 
be calculated. Figure ([1]) shows the entanglement of two spins, one at site 1 and the 
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30 0.02 



Figure 1: (Color Online) The entanglement of two spins at sites 1 and r as a function 
of their distance and the coupling g. Near the critical point the range of entanglement 
increases at the cost of its amplitude, for every g there is a maximum distant, beyond 
which there is no entanglement. 



other at site r as a function of r and the parameter g. It is seen that in these models the 
range of entanglement increases as we approach the critical point, however its value 
decreases when we approach this point. We can obtain an approximate expression 
for the maximum range of entanglement. Inspection shows that the only eigenvalue 
of p'^^{l,r) which may go negative is 



+ 7 1 F/ v5 7r~^ 
= ^ -^J [a - -fY + Sfi'' , 

note that = in the thermodynamic limit (llSp . Thus the inequality < 
determines the range of entanglement. This is however equivalent to 07 < 2fj? or in 
view of (fT8]) 

{l + 2\g\Y-'-{l-2\g\Y-'<2{l + 2\g\f-^. 

For small values of g, i.e. \g\ <C 1, we can write 1 it 2\g\ ~ e^^'^' and then the above 
inequality transforms to 

/n(3) 

implying that a entanglement exists up to a range of vq ~ + 2. As an example 
for g = 0.02, this gives tq ~ 15, which is also evident from figure ([1]). 
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Figure 2: The single site entanglement for models with a = 1 (left) and a = — 1 
(right), for several values of system sizes, N=15, 21, 25, 31, and 35. For a = 1, the 
entanglement is almost independent of size, as long as > 10 and for cr = —1, it is 
so unless g is very close to the critical point. 

5 Systems with finite size 

The interest in entanglement properties of spin systems stems not only from its pos- 
sible relation to the critical properties of such systems, but also from their possible 
candidacy for the future implementation of quantum computers. In this case we are 
dealing with a finite array of interacting spins which has relaxed to its ground state. 
Tuning the interactions of these spins with each other, changes the ground state and 
it is desirable to have controllable entanglement between different spins of this array. 
Therefore in this section we study more closely the properties of such matrix product 
states for finite values of A^. 



5.1 Single site entanglement 

To measure the entanglement of one site with the rest of the chain we calculate the 
von Neumann entropy of the density matrix of one site, which is readily found to be 




(20) 



with 



a = 




6 = 1 - 2a 



(21) 



Z 



and 




(22) 
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Figure 3: (Color Online) The entanglement of two distant spins as a function of their 
distance for a ring of size = 40 for several values of g near the critical point. 

The one-site entropy which measures the entanglement of one site with the rest of 
the lattice will then be given by 



For >> 1 one can verify that a — 5, therefore for system sizes, S{g, N) is indepen- 
dent of N. Figure ([2]) shows the behavior of one-site entanglement for different system 
sizes. It shows clearly the difference between the models with a = 1 and a = —1. 
Note that the difference of cr = —1 and a = 1 models shows up only for systems of 
odd size. 

5.2 Entanglement of two distant sites 

The entanglement of two sites can be measured by the negativity of the reduced den- 
sity matrix of two sites. The general form of p(l, r) is the same as in (I17p except that 
the correlation functions now depend on the system size in a rather complicated 
way. We do not write the explicit form of for finite and only report the 

behavior of entanglement in figures ([3]) and Figure ([3]) shows the entanglement 
of two distant spins as a function of their distance for a ring of size A^ = 40 for several 
values of g near the critical point. Several features are evident. First the entangle- 
ment has always a finite range. Furthermore, the range of entanglement increases as 
we approach the critical point, however its value decreases. Near the critical point, it 
is well known that correlation lengths diverge, here we see that entanglement range 
also diverges although this is accompanied by lowering of its value. Figure (jH) shows 
entanglement between adjacent spins as a function of the size of the system at a fixed 
value of g. The figure shows that two spins at distance 4 (r = 5) can be entangled 
for rings of size up to A^ = 26. Also for any value of g and any distance r, there is 
a maximum system size Nmax{f,g) above which those two sites can not be entangled 
at all. This figure shows that Nmax{f,g) decreases with r. 



S 



2alog(a) - (1 - 2a)log(l - 2a). 



(23) 



9 



Negativity 




f 10 20 30 40 50 60 

N 



Figure 4: (Color Online) The entanglement of two spins with distances 1,2,3 and 4 
(from top to bottom) as a function of the system size. 

In figure ([5]) the entanglement of nearest neighbor sites is displayed for several 
values of system sizes when a = 1. At the critical point, the nearest neighbor spins 
are not entangled, however, entanglement appears for any infinitesimal deviation from 
this point. On the other hand when a = —1 (the class of models which contain the 
AKLT point), there is entanglement at g = 0, only for rings of odd size (figure [6]). 
The behavior for even-sized rings is identical with the case when a = 1. 

Figure ([5]) suggests a scaling behavior for entanglement near the critical point. To 
investigate this property, we consider for definiteness, the case a = 1. For each value 
of N, the entanglement attains a maximum at a point gm{N), where its value at this 
point is denoted by £m{N)- In the inset of figure ([7]) we plot log{gm) and log{£m) 
versus log (N). We find numerically that 

log(c,^) = -1.0771og(iV)- 0.106 Ri- log A^-c, 

log(£:„) = -1.0861ogA-0.214Ri -logA-d, (24) 

where c and d are two constants, independent of N. From the scaling hypothesis, this 
means that one can write the negativity in the vicinity of the critical point as 

£{g,N)^j^f{Ng), 

where / is a universal function. In figure ([7j) we plot log(A<S) as a function of log{Ng). 
It shows that all the data collapse on a single curve for 4 order of magnitudes of gN. 

5.3 Limiting form of the states 

It is instructive to find the explicit form of the state near the critical point, i.e. when 
\g\ <C 1. For such an analysis we should find the dominant amplitudes in the 
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Figure 5: (Color Online) The entanglement of adjacent spins as measured by their 
negativity for different values of system size, from top to bottom equal to A'^ = 
15,20,25,30,35,00(0- = 1). 



linear superposition of all states. We consider the cases a = 1 and a = —1 separately. 

Case a: a=l Near the critical point \g\ <C 1, the dominant amplitudes are ipoo-.-o, 
ipf^j and ip-j^j^, where ip^^j denotes the amplitude of a state \k,l) in which two spins at 
sites k and I are respectively excited to 1 and —1. To first order in g, the state ([T]) 
becomes 

N 

\^)c^\0,Or--0)-g ^ i\k,l) + \k,l)). (25) 

k<l=l 

It is not difficult to find the negativity of this state which is 

£{p{l,r))=2\g\. (26) 

Case b: cj=-l, N=even For \g\ <^ 1, the dominant amplitudes are V'oo - o oc 1 and 

^|J,J oc tr{A'^-'AiA',~'^''AjA^-') « tr{A,A',-''~'A-) = 25(-l)'-^ 

and ■0^; oc 25f(— 1)'^"', where we have used the fact that ^q^vIi = Ai for any m. 
Thus the state becomes 

1^) :^ |0, 0, • • • 0) + gJ2{-iy-\\k:i) + \k, I)). (27) 

k<l 

From the two-site density matrix, one can find the negativity which is given by 
£{p{l,r)) = 2\g\. 

Case c: a=-l, N=odd In this case we have V^oo- o = and so the state becomes 
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Figure 6: (Color Online) The entanglement of adjacent spins as measured by their 
negativity for different systems with odd-size, from top to bottom equal to = 
15,21,25,31,35, {a = -1). 



The entanglement of such a state can be calculated by determining the reduced 
density matrix of two sites. After a rather lengthy but calculation, we find 

W,r)) = ^^v(F3I)l(^-2)(^-3)-l 



[{N - 2){N - 3) + 1]2 + 8{N - 2)2| 



6 Conclusion 

In this paper we have studied in detail, the entanglement properties of a general 
one-parameter family of matrix product state of spin one particles defined on a ring. 
The state has some plausible symmetries, like rotational symmetry around the z, 
axis, and the parity symmetry. Such a state is the ground state of local Hamiltonian 
defining the nearest neighbor interaction of spins. The state goes a sharp transition 
when its continuous parameter, denoted by g, passes through a critical point. In the 
thermodynamic limit, this can be ascribed to a quantum phase transition of the system 
described by the local Hamiltonian. We have studied the entanglement properties of 
the state, near this transition point, both in the thermodynamic limit and for finite 
chains. The study of finite chains is motivated by the possible role of such systems in 
quantum computing. We have considered two measures of entanglement, namely the 
entanglement of one site with the rest of the chain, (also studied in [T8] for infinite 
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Figure 7: (Color Online) Near the critical point, the entanglement of nearest neighbor 
sites, as measured by their negativity, shows a scaling behavior. The main plot shows 
log{N£) versus log{gN), for values of iV = 20, 25, 30 • • • 80, 85. Ah the data collapse 
on one single curve, for a range of 4 orders of magnitude in the value of gN. 

rings) and the entanglement of two distant sites with each other. Our findings can be 
summarized as follows: 1- When the parameter g approaches its critical value go = 0, 
the range of entanglement (between distant spins) increases indefinitely at the cost 
of its magnitude, fig. ([T]), 2- for any non-zero value of g, there is a threshold distance 
beyond which there is no entanglement between spins, fig. ([3]), 3- for any two spins 
with a fixed distance, there is a threshold system size, above which entanglement 
vanishes, fig. (|4|), and finally there is a kind of scaling behavior in entanglement 
properties of neighboring sites, near the critical point (Figure [Tj). 
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